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This paper describes an efficient and unique method for computing the shortwave di- 
rect radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using 
CALIOP and MODIS data. It accounts for the overlapping of aerosol and cloud rigor- 
s ously by utilizing the joint histogram of cloud optical depth and cloud top pressure. 
Effects of sub-grid scale cloud and aerosol variations on DRE are accounted for. It is 
computationally efficient through using grid-level cloud and aerosol statistics, instead of 
pixel-level products, and a pre-computed look-up table in radiative transfer calculations. 
We verified that for smoke over the southeast Atlantic Ocean the method yields a sea- 
io sonal mean instantaneous shortwave DRE that generally agrees with more rigorous 
pixel-level computation within 4%. We have also computed the annual mean instanta- 
neous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over 
global ocean based on 4yr of CALIOP and MODIS data. We found that the variability 
of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly 
is driven by the optical depth of the underlying clouds. 

1 Introduction 

The shortwave direct radiative effect (DRE) of aerosols at the top of the atmosphere 
(TOA) is strongly dependent on the reflectance of the underlying surface. Over dark 
surfaces (e.g. ocean, vegetated land), the scattering effect of aerosols is generally 
20 dominant, leading to negative DRE (i.e., cooling) at TOA (Yu et al., 2006). In con- 
trast, when light-absorbing aerosols occur above clouds or other bright surfaces (such 
as snow, ice, desert), aerosol absorption is significantly amplified by cloud or surface 
reflection, offsetting or even exceeding the scattering effect of the aerosol leading to 
a less negative or even positive (i.e., warming) TOA DRE (Abel et al., 2005; Keil and 
25 Haywood, 2003; Twomey, 1977). Therefore, in order to understand the full complex- 
ity of aerosol radiative effects on climate, it is important to quantify the DRE under 
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both clear-sky and cloudy-sky conditions. Although the DRE of aerosols in clear-sky 
regions has been extensively studied and is relatively well constrained based on ad- 
vanced satellite remote sensing measurements acquired in the last decade (Yu et al., 
2006), the cloudy-sky DRE is generally assumed to be negligible or simulated by mod- 
5 els (Schulz et al., 2006). Current model simulations show a large inter-model spread in 
cloudy-sky DRE (Schulz et al., 2006), which results from inter-model differences in both 
aerosol and cloud properties (Schulz et al., 2006; Stier et al., 2013). There is a clear 
need for an observational constraint on the DRE of above-cloud aerosol (ACA). 

Recent advances in satellite remote sensing techniques have provided an unprece- 
io dented opportunity for studying the DRE of ACA. In particular, the availability of mea- 
surements from the space-borne Cloud-Aerosol Lidar with Orthogonal Polarization 
(CALIOP) sensor onboard NASA’s Cloud-Aerosol Lidar and Infrared Pathfinder Satel- 
lite Observations (CALIPSO) satellite has provided a revolutionary global view of the 
vertical distribution of aerosols and clouds (e.g., Winker et al., 2013). Using CALIOP 
is aerosol and cloud layer products, Devasthale and Thomas (2011) found frequent oc- 
currences of aerosols residing above low-level clouds in several regions of the globe. 
In particular, they found a high frequency of smoke occurrence over low clouds in the 
southeast Atlantic, western coasts of South America (e.g., Columbia, Ecuador, and 
Peru) and southern Asia. These authors also found a high frequency of natural and 
20 polluted dust aerosols overlapping low clouds off the western coasts of Saharan Africa 
in boreal summer and over the eastern coast of China in boreal spring (see Fig. 3 of 
Devasthale and Thomas, 2011). 

CALIOP measurements of ACA properties, in combination with satellite cloud prod- 
ucts from, for example, the Moderate Resolution Imaging Spectroradiometer (MODIS), 
25 have been used in several recent studies to derive the DRE of ACA with radiative 
transfer simulations (e.g., Chand et al., 2009; Costantino and Breon, 2013b; Oikawa 
et al., 2013). Chand et al. (2009) used CALIOP above-cloud AOD retrievals (Chand 
et al., 2008) and Terra-MODIS cloud products, both aggregated to 5° gridded monthly 
means, to calculate the radiative effects of smoke transported above the low-level 
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stratocumulus deck in the southeastern Atlantic. A major point made in this study is 
that the all-sky DRE of elevated light-absorbing aerosols, such as transported smoke, 
is strongly modulated by the underlying cloud properties. However, the spatial-temporal 
aggregation of both CALIOP and MODIS data to coarse gridded monthly means ob- 
5 scures the potential influence of cloud and aerosol variability on the DRE. In particular, 
using grid box mean cloud optical depth for DRE calculation might lead to biases in 
DRE due to the plane-parallel albedo bias (Oreopoulos et al., 2007; Min and Zhang, 
2013). Moreover, the MODIS level-3 aggregation algorithm samples all liquid water 
clouds, regardless of possible retrieval contamination by ACA. As a result, the total 
io population of liquid water clouds in the MODIS level-3 products (daily or monthly) may 
be significantly different from that of below-aerosol-only cloud population. Therefore, 
using level-3 MODIS products without distinguishing below-aerosol-only cloud popu- 
lation from the total can potentially lead to significant errors. The problem could be 
further complicated by biases in MODIS cloud retrievals associated with the presence 
is of overlying light-absorbing aerosols. When a cloud-pixel is contaminated by overlying 
light-absorbing aerosols the MODIS cloud optical depth (COD) retrieval is generally 
biased low (Coddington et al., 2010; Haywood et al., 2004), an effect not considered in 
most previous studies (e.g., Chand et al., 2009; Costantino and Breon, 2013b; Oikawa 
et al., 2013). Meyer et al. (2013) found that correcting the MODIS COD bias due to 
20 ACA contamination can lead to a more positive ACA DRE. 

Most recently, Meyer et al. (2013) collocated CALIOP above-cloud AOD and Aqua- 
MODIS cloud properties at pixel level, and the DRE was then computed at individual 
collocated pixels. Such rigorous collocation has obvious advantages as it takes into 
account the sub-grid variability of clouds and aerosols. However, such a method is 
25 computationally expensive and requires large amounts of pixel-level data, which poses 
a major challenge for application on a global scale and multiyear basis. 

The objective of this paper is to describe a novel method for computing the DRE 
of ACA. In this method, we attempt to balance the need for computational efficiency 
with the need for rigorous treatment of aerosol-cloud overlapping and small-scale 
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variabilities of aerosol and clouds. Our method has several unique features: (1) it takes 
sub-grid scale cloud and aerosol variation into account in DRE computation; (2) it treats 
the overlapping of aerosol and cloud rigorously by utilizing the joint histogram of cloud 
optical depth and cloud top pressure in the MODIS level-3 product; (3) it is computa- 
s tionally efficient because of the use of a pre-computed look-up table of ACA DRE. In the 
following sections, we briefly introduce the CALIOP and MODIS data used (Sect. 2), 
describe the key assumption and features of the method (Sect. 3), validate it through 
comparison with pixel-level computations as in Meyer et al. (2013) (Sect. 4), and con- 
clude with a summary and discussion (Sect. 5). 


io 2 Satellite data 

2.1 CALIOP aerosol and cloud layer products 

Since its launch in 2006, the space-borne lidar CALIOP has continuously acquired, with 
near global (albeit instantaneously sparse) coverage, attenuated backscatter measure- 
ments at 532 nm and 1064nm, including linear depolarization information at 532 nm 
is (Winker et al., 2009). The CALIOP level-2 retrieval algorithm consists of several steps. 
First, a “feature finder” algorithm and cloud-aerosol discrimination (CAD) algorithm are 
used to detect aerosol and cloud layers, and record their top and bottom heights and 
layer integrated properties (Vaughan et al., 2009). Second, the detected aerosol lay- 
ers are further classified into six sub-types (i.e., polluted continental, biomass burning, 
20 desert dust, polluted dust, clean continental and marine) (Omar et al., 2009) and cloud 
layers into different thermodynamic phases (Hu et al., 2007) based on the observed 
backscatter, color ratio and depolarization ratio. Third, a priori lidar ratios, selected 
based on aerosol sub-type and cloud phase, are used to derive the extinction of an 
aerosol or cloud layer from the attenuated backscatter profile (Young and Vaughan, 

25 2008). 
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In this study, we use CALIOP level-2 version 3 aerosol and cloud layer prod- 
ucts at a nominal 5 km horizontal resolution (i.e., CAL_LID_L2_05kmALay and 
CAL_LID_L2_05kmCLay) for aerosol-cloud overlapping detection, and for information 
on aerosol layer properties, including aerosol layer type, aerosol optical depth (AOD), 
s and layer top and bottom height. In addition to physical properties, the CALIOP layer 
products also provide various metrics and flags for data quality assurance. These in- 
clude CAD score (Liu et al., 2009), horizontal averaging scale, the extinction quality 
control (QC) flag, and estimated uncertainty for layer AOD. In this study, we use these 
metrics following the best practice advice of the CALIPSO science team to screen for 
io reliable retrievals (e.g., Winker et al., 2013) (see Table 1). 

2.2 MODIS daily level-3 cloud property product 

We use the Collection 5 (C5) Aqua MODIS level-3 daily gridded Atmosphere product 
(i.e., MYD08_D3) for cloud properties and other parameters, such as solar zenith an- 
gle, needed for ACA DRE computations. The MODIS level-3 gridded products are ag- 
15 gregated at 1° resolution from the MODIS level-2 pixel-level retrievals (Hubanks et al., 
2008; King et al., 2003). As summarized in Platnick et al. (2003), the operational level- 
2 MODIS cloud product provides cloud masking (Ackerman et al., 1998), cloud top 
height retrieval based on C0 2 slicing or the infrared window method (Menzel et al., 
1983), cloud top thermodynamic phase determination (Menzel et al., 2006), and cloud 
20 optical and microphysical property retrieval based on the bi-spectral solar reflectance 
method (Nakajima and King, 1990). In addition to these cloud parameters, the level-2 
products also provide pixel-level runtime Quality Assessment (QA) information, which 
includes product quality as well as processing path information. 

All MODIS level-2 atmosphere products, including the cloud, aerosol and water vapor 
25 products, are aggregated to 1° spatial resolution on a daily, eight-day, and monthly 
basis. Aggregations include a variety of scalar statistical information (mean, standard 
deviation, max/min occurrences) and histograms (marginal and joint). A particularly 
useful level-3 cloud product for this study is the daily joint histogram of cloud optical 
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depth (COD) vs. cloud top pressure (CTP), derived using daily counts of successful 
level-2 pixel retrievals that fall into each joint COD-CTP bin. Eleven COD bins, ranging 
from 0 to 1 00, and 1 3 CTP bins, ranging from 200 to 1 000 mb, comprise the histogram. 
As discussed below, the COD-CTP joint histogram allows for identifying the portion of 
5 the cloud population that lies beneath the aerosol layer found by CALIOP, as well as the 
corresponding COD probability distribution needed for DRE estimation. In addition to 
the COD-CTP joint histogram, we also use the gridded mean solar and sensor zenith 
angles for calculating DRE and correcting COD bias associated with ACA. 


3 Methodology 
io 3.1 Theoretical basis 

As in previous investigations (e.g., Chand et al., 2008, 2009; Costantino and Breon, 
2013b; Meyer et al., 2013), we focus on the simplest case of overlapping aerosol and 
cloud, i.e., a single layer of aerosol overlying a single layer of low-level liquid-phase 
clouds, which is commonly observed in many regions of the globe (Devasthale and 
is Thomas, 2011). More complex situations exist, such as an aerosol layer located in be- 
tween high cloud and low cloud, or an aerosol layer overlying multiple layers of clouds. 
However, such situations are either beyond the capability of CALIOP or relatively rare 
(Devasthale and Thomas, 2011). As such, they are left for future research and are not 
considered here. 

20 To illustrate the theoretical foundation of the method, consider the schematic ex- 
ample in Fig. 1. For a given grid box, the gridded mean broadband shortwave DRE 
((DRE) aca ) averaged over all ACA pixels within the grid box is given by: 

CO CO 

(DRE) aca = JJ DRE(T c ,T a )p(T c ,T a )dT c dT a , (1) 


where p(r c ,r a ) is the joint probability density function (PDF) of the above-cloud AOD at 
532 nm (r a ) and below-aerosol COD (t c ) of ACA pixels. We note that, in addition to r a , 
DRE is also a function of wavelength dependence of AOD, spectral single scattering 
albedo and asymmetry factor, which are not explicitly included in this equation. These 
s properties are computed using a Mie scattering code (Wiscombe, 1980) based on the 
aerosol model described in Meyer et al. (2013). The dependencies on solar zenith 
angle, surface reflectance, cloud particle effective radius, and atmospheric profile are 
also omitted; solar zenith angle and surface reflectance are expected to have only 
minor variation within the grid box, while the impact of cloud particle effective radius 
io and atmospheric profile on shortwave DRE is relatively small. Since p(r c , r a ) describes 
the covariation of aerosols and clouds for the ACA pixels, it should ideally be derived 
from collocated CALIOP aerosol and MODIS cloud retrievals at pixel level as in Meyer 
et al. (2013). This requires large amounts of pixel-level data, however, as one month 
of global daytime C5 MODIS level-2 cloud products in HDF format are roughly 150 Gb. 
is Therefore, pixel-level collocation and radiative transfer simulation are too computation- 
ally expensive and cumbersome for multiyear global studies. 

A key assumption in our method, which allows us to avoid tedious pixel-level collo- 
cation, is that above-cloud AOD and below-aerosol COD are statistically independent. 
Under this assumption, p( T c ,r a ) = p(r c )-p( r a ) and Eq. (1) reduces to: 


20 


(DRE) ACA 
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P(T a )dT a 
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where p(r c ) and p(r a ) are the PDF of below-aerosol COD t c and above-cloud AOD 
T a , respectively, of ACA pixels. The advantage of Eq. (2) is that it allows p(r c ) and 
p(r a ) to be derived separately and independently. This assumption is reasonable con- 
sidering that transported ACAs and low-level boundary layer clouds are usually well 
25 separated vertically (Devasthale and Thomas, 201 1 ) and controlled by different meteo- 
rological conditions. The potential coupling between the two is that overlying absorbing 

10000 


Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper 


aerosols could influence the evolution of clouds through changing atmospheric strati- 
fication (Wilcox, 2010). However, a recent observational study (Costantino and Breon, 
2013a) found no correlation between above-cloud AOD and below-aerosol COD. More- 
over, as shown in Sect. 4, we have compared the DRE derived from pixel-level colloca- 
5 tion (i.e., based on Eq. 1) with that from independent sampling of p(r c ) and p( r a ) (i.e., 
based on Eq. 2) and found very good agreement. 

In our method, the PDF of above-cloud AOD p(r a ) is derived from the CALIOP 5 km 
aerosol and cloud layer products through the following steps: (1 ) for each 5 km CALIOP 
profile that falls within a given latitude-longitude grid box, we first search for an aerosol 
io layer; ( 2 ) if an aerosol layer is detected and the quality metrics pass the quality as- 
surance criteria summarized in Table 1, we then proceed to check for the presence of 
an underlying liquid-phase cloud layer within the profile using the CALIOP cloud layer 
product; (3) if a cloud layer is present, the AOD of the aerosol layer is recorded for the 
derivation of the p(z a ) of the grid box. The bottom height of the aerosol layer is also 
is recorded to derive the grid mean aerosol layer bottom height. Once all of the CALIOP 
profiles within the grid box are processed, we obtain the PDF of the above-cloud AOD 
p(r a ) and the mean aerosol layer bottom pressure (P bottom ). 

As schematically illustrated in Fig. 1 , the PDF of below-aerosol COD p(r c ) is derived 
from the joint histogram of cloud optical depth and cloud top pressure (COD-CTP joint 
20 histogram) in the MODIS daily level-3 product, using the grid mean aerosol layer bottom 
pressure (p b0tt0 m) derived above. For a given grid box, we first identify the population 
of liquid-phase clouds below the pressure level (bottom)- This subset, together with the 
AOD PDF p(r a ), is then used to calculate DRE according to Eq. (2). 

3.2 DRE look-up tables 

25 For better computational efficiency, we use pre-computed aerosol-type specified look- 
up-tables (LUTs), instead of online radiative transfer computation, when deriving the 
DRE of ACA. The concept of our LUTs is somewhat similar to the “radiative kernels” 
described in Hartmann et al. (2001) and Zelinka et al. (2012) for computing cloud 
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radiative feedbacks. The LUT for each aerosol type consists of DREs at both TOA and 
surface (not used in this study) for various combinations of AOD, COD, CTP and solar 
zenith conditions. As such, once the aerosol type and AOD are known from CALIOP 
and COD, CTP and solar zenith angle are known from MODIS, the corresponding 
s DRE can be obtained through LUT interpolation. Note that the CALIOP only provides 
AOD at lidar wavelengths (e.g., 532 nm and 1064nm) for each aerosol type. Therefore, 
broadband aerosol scattering properties, including spectrally dependent AOD, single- 
scattering albedo and asymmetry factor, are needed for the development of LUT. The 
current version of LUT focuses on light-absorbing aerosols (e.g., smoke and polluted 
io dust). In order to validate our method with more rigorous pixel-level computations, we 
adopt the broadband aerosol optical properties from Meyer et al. (2013) in the com- 
putation of the current LUT. The aerosol model in Meyer et al. (2013) is extended 
from an absorbing aerosol model developed for the MODIS Collection 5 Aerosol Prod- 
uct (MOD04) (see Table 4 of Levy et al., 2009). The MOD04 aerosol models define 
is aerosol size distributions and refractive indices dependent solely on prescribed AOD 
at 550 nm (MODIS band 4; note that the absorbing aerosol model used here assumes 
a constant index of refraction, 1 .51 -0.02/', at all wavelengths). At AOD = 0.5 at 550 nm, 
the single-scattering albedo of this model is about 0.9 over the visible spectral region 
(see Fig. 7 of Meyer et al., 2013), which is in the range of previously reported values 
20 (e.g., Keil and Haywood, 2003; Myhre et al., 2003). The current AOD bins (at 550 nm) in 
the LUT range from 0.05 to 1 .5, which covers most of the above-cloud AOD observed 
by CALIOP. The current COD bins, logarithmically spaced, range from 0.1 to 300. Fol- 
lowing the MODIS level-3 data, the thirteen CTP bins range from 1000 mb to 200 mb. 
The solar zenith angle bins range from 0 to 80°. Radiative transfer computations are 
25 carried out using the RRTM-SW model (Clough et al., 2005; lacono et al., 2008). Lam- 
bertian ocean surface reflectance is set to 5%. Cloud droplet effective radius is fixed 
at 15 pm, which is close to the global mean value over oceans observed by MODIS 
(King et al., 2013). Water cloud optical properties are calculated internally by RRTM. 
For atmospheric profiles of water vapor and temperature, we use NCEP R1 reanalysis 
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data (Kistler et al., 2001) averaged both zonally and annually. Our sensitivity tests in- 
dicate that the shortwave DRE of ACA is very insensitive to cloud effective radius or 
atmospheric profiles. 

3.3 Cloud optical depth correction 

5 As noted in previous studies (Coddington et al., 2010; e.g., Haywood et al., 2004), when 
a cloudy MODIS pixel is contaminated by overlying light-absorbing aerosols the COD 
retrieval is generally biased low. We have developed a fast COD correction scheme to 
account for the COD retrieval bias due to ACA in our DRE computation, which is illus- 
trated in Fig. 2. This scheme requires both the cloud reflectance LUT for clouds without 
io ACA, for which we use the MODIS operational LUT, and clouds with ACA, for which we 
use the one developed by Meyer et al. (2013). In the operational MODIS retrieval, the 
reflectance LUT of cloud without ACA is used to interpret the reflectance of all clouds, 
including those affected by ACA. Based on this fact, we first infer the “observed” cloud 
reflectance (after atmospheric correction) by interpolating the reflectance LUT of cloud 
is without ACA corresponding to the biased COD. Then, we use the “observed” cloud 
reflectance and ACA-affected LUT (derived based on CALIOP AOD) to determine the 
corrected COD. This COD correction process is performed for every combination of 
COD bin in p(r c ) and AOD bin in p(r a ). In the final step we resample the corrected 
CODs to obtain the corrected p(r c ). 

20 It should be noted that because different aerosol types may have different impacts 
on MODIS COD retrievals, the above COD correction process is aerosol-type depen- 
dent. In this study, we use light-absorbing aerosols as an example to illustrate our 
method, and for validation purposes we use the aerosol model developed in Meyer 
et al. (2013) for the development of LUTs for DRE computation and COD correction. 
25 However, the LUTs can be easily extended to other aerosol models. In fact, as part of 
ongoing research, we are extending our LUTs to include all six operational CALIOP 
aerosol models as described in (Omar et al., 2009). 
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4 Implementation and validation of new DRE estimation scheme 

Each year during austral winter, dry season biomass burning activities throughout 
southern Africa inject large amounts of smoke into the troposphere (Eck et al., 2003; 
Ichoku et al., 2003; Myhre et al., 2003). Prevailing easterly winds during this season 
s often transport the smoke westward off the continent, over the ocean, where exten- 
sive marine boundary layer clouds persist for most of the year. Under the descending 
branch of the Hadley cell, the air mass above the boundary layer is quite dry. Due to the 
lack of efficient wet scavenging, the transported aerosol layers can remain suspended 
in the atmosphere for days, creating a near-persistent smoke layer above the stratocu- 
io mulus deck over the southeastern Atlantic Ocean during this time (Chand et al., 2009; 
Devasthale and Thomas, 2011; Keil and Haywood, 2003; Wilcox, 2010). 

To validate our method, we have compared the DRE of above-cloud light-absorbing 
aerosols in this region with pixel-level computations from Meyer et al. (201 3). Figure 3a 
shows the seasonal mean (August/September 2007-2011) instantaneous TOA DRE 
is of above-cloud smoke and polluted dust based on the pixel-level computations from 
Meyer et al. (2013). Figure 3b shows the corresponding instantaneous TOA aerosol 
radiative forcing efficiency (RFE) defined as the DRE per unit AOD. The DRE and RFE 
results computed using our method described in the previous section are shown in 
Fig. 3c and d, respectively. Evidently, both DRE and RFE computed using our new 
20 method agree closely with the pixel-level computations. Figure 4 shows the meridional 
mean DRE and RFE for the region based on the results in Fig. 3. Not surprisingly, the 
results based on the two methods are almost identical. Note that the CODs used in the 
computations for Figs. 3 and 4 are directly from the MODIS products without COD cor- 
rection. We have also compared the DRE and RFE from the two methods based on the 
25 corrected COD and achieved again very good agreement (these results are not shown 
here because they are quite similar to those in Figs. 3 and 4). The seasonal and re- 
gional mean DRE and RFE, based on the corrected COD, from the pixel-level computa- 
tion method in Meyer et al. (201 3) are 6.63 W m -2 and 55.97 W rrT 2 AOD -1 , respectively 
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(see Table 2). The corresponding values from our new method are 6.39 Wm " 2 and 
53.77 WirT 2 AOD -1 , respectively. As shown clearly in Figs. 3 and 4 and Table 2, the 
DRE inferred from our new method agrees very well with the pixel-level computations. 
Furthermore, the difference between the two methods is much smaller than, for ex- 
5 ample, the uncertainty associated with CALIOP retrieval biases. Recent studies found 
about a factor of two difference between CALIOP nighttime and daytime AOD retrievals 
(Chand et al., 2008; Meyer etal., 2013). If this difference is considered as CALIOP AOD 
retrieval uncertainty, it could lead to 30 %~ 50% uncertainty in ACA DRE computation 
(see Fig. 12 of Meyer et al., 2013). 


io 5 Summary and discussion 

Recent advances in satellite-based remote sensing, in particular the launch of the 
space-borne lidar CALIOP, have provided an unprecedented opportunity for studying 
the radiative effects of above-cloud aerosol (ACA). However, the methodologies used 
in recent studies for computing the ACA DRE appear to be either oversimplified (e.g., 
is Chand et al., 2009; Oikawa et al., 2013) or too cumbersome (e.g., Meyer et al., 2013). 
This paper describes a novel method recently developed for computing the shortwave 
DRE of above-cloud aerosols over ocean. Our method has several unique features 
compared to methods used in previous studies: ( 1 ) it takes sub-grid scale cloud and 
aerosol variation into account in DRE computations, similar to Meyer et al. (2013); (2) 
20 it also treats the overlapping of aerosol and cloud rigorously by utilizing the joint his- 
togram of COD and CTP in the MODIS level-3 cloud product; (3) it differs from Meyer 
et al. (2013) in its reliance on grid-level cloud statistics (i.e., COD-CTP joint histogram), 
instead of pixel-level products, and utilizes pre-computed look-up tables for ACA DRE 
computations, making it thus much more efficient than pixel-level computations. As 
25 shown in Figs. 3 and 4 and Table 2, DRE computed using our method agrees well with 
the pixel-level computations in Meyer et al. (2013). 
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In addition to the Southeast Atlantic region, we have recently begun investigating 
the DRE of above-cloud light-absorbing aerosols for global ocean. Some preliminary 
results are shown in Fig. 5. We first derived the daily grid-level statistics of above-cloud 
AOD and below cloud COD, as well as the corresponding ACA DRE, using the method 
s described above and then aggregated the daily means to annual mean. The temporal 
aggregation is weighted by the number of ACA pixels in each day during 2007-2010. 
For example, the annual mean ACA DRE in Fig. 5c is aggregated from daily mean 
based on the following equation: 

Ia/,-<dre,) aca 

( dre ) aca = - — jjj: . (3) 

/ 

io where (DREj ) aca is mean instantaneous ACA DRE in each day averaged over ACA 
pixels, A/ ( is the number of ACA pixels in the grid box in each day, and (DRE) ACA is 
the annual mean instantaneous ACA DRE shown in Fig. 5c. In Fig. 5a a global map is 
shown of the annual mean 550 nm AOD of above cloud smoke and polluted dust de- 
rived based on 4yr (2007-2010) of CALIOP aerosol and cloud layer products. Similar 
is to Devasthale and Thomas (201 1 ), we note several “hotspots” of ACA over the South- 
east Atlantic, the East-Central Atlantic off the western coast of Saharan Africa, the 
Arabian sea, and the North Pacific basin off the coast of eastern Asia. It is interesting 
to note that the ACA AOD over the east-central Atlantic and Arabian Sea is noticeably 
larger than that over the southeast Atlantic and North Pacific basin. Figure 5b shows 
20 the annual mean below-aerosol COD derived from the MODIS daily level-3 cloud prod- 
uct using the method described in Sect. 3. A notable feature in the figure is that the 
below-aerosol COD over the North Pacific basin is significantly larger than that over 
other ACA regions. Figure 5c shows the annual mean shortwave DRE at TOA aggre- 
gated from daily values due to ACA smoke and polluted dust over the global ocean. 
25 It is intriguing to see that the DRE of ACA over the North Pacific basin is significantly 
larger than that over the southeast Atlantic, which is in turn larger than the DRE over 
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the east-central Atlantic and the Arabian Sea. In fact, some negative DREs are ob- 
served in the latter two regions. This is interesting because the above-cloud AOD over 
these regions is actually larger, while the below-aerosol COD over these regions is 
smaller, in comparison with the counterparts over the southeast Atlantic and North Pa- 
5 cific basin. Therefore, the preliminary results seem to suggest that the variability of 
DRE of ACA is modulated by COD, rather than AOD, although it should be noted that 
we have focused only on the light-absorbing aerosols, i.e., smoke and polluted dust, 
and assumed the same broadband scattering properties for them as in Meyer et al. 
(2013). Future research is needed to study the impact of aerosol type and scattering 
10 properties on the temporal-spatial variation of DRE on a global scale. Nevertheless, 
the preliminary results shown in Fig. 5 clearly demonstrate the usefulness of our new 
method for global study. 

It should be noted that this study, and previous ones using CALIOP observations 
(e.g., Chand et al., 2008; Meyer et al., 2013; Oikawa et al., 2013), are limited by the ca- 
is pabilities of CALIOP. Arguably, some aerosols exist above every cloud. However, not all 
ACA can be detected by CALIOP due to its inherent limitations. Some ACAs are sim- 
ply too optically thin to be detected, though their radiative effects are also expected to 
be small. There are also other possibilities. For example, a confined aerosol layer has 
larger volume backscatter than a vertically stretched layer, even if the aerosol amounts 
20 are the same, and therefore is more easily detected by CALIOP. Passive sensors, on 
the other hand, are less affected by the vertical distribution of ACA because they ob- 
serve column-integrated scattering by aerosols. Recently, several novel techniques 
have been developed to detect and retrieve ACA properties using passive sensors. 
Waquet et al. (2009) developed a method based on multi-angular polarization mea- 
25 surements from POLDER (Polarization and Directionality of the Earth Reflectances) to 
retrieve the AOD of above-cloud smoke. This method has recently been extended to 
include both smoke and dust aerosols (Waquet et al., 201 3). Most recently, Jethva et al. 
(2013) demonstrated a color ratio method to retrieve the above-cloud AOD based on 
MODIS multiple spectral cloud reflectance measurements. A review of the emerging 
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satellite-based observations of above-cloud aerosols can be found in Yu and Zhang 
(2013). The capabilities and limitations of the passive techniques need to be system- 
atically studied through inter-comparison and comparison with CALIOP observations, 
but they may provide a complementary perspective on ACA. Note that passive imagers 
s have much larger spatial coverage than CALIOP, which makes calculations of the DRE 
at the pixel level computationally expensive. In this regard, our method satisfies the 
need for efficiency of ACA DRE computations based on passive imager retrievals. 

As a final remark, we would like to point out that the ACA DRE discussed in this 
study is still a few steps away from the all-sky aerosol radiative effect ((DRE) a|| . sky ). For 
io a given grid box, the (DRE) a|| _ sky can be decomposed into the sum of clear-sky and 
cloudy-sky DRE: 

<DRE) al , sky = (1 - g • (DRE) c|ear + f c . f ACA . (DRE) aca , (4) 

where f c is the cloud fraction, (DRE) clear is the DRE averaged over the clear-sky portion 
of the grid box, f ACA is the fraction of cloudy pixels with ACA detected by CALIOP or 
is other sensors, and (DRE) ACA is the DRE averaged over all ACA containing pixels. It 
is important to note an implicit assumption made in Eq. (4), that is, when a distinct 
ACA layer is not detected, the DRE is zero. Different sensors (or different retrieval 
algorithms for the same sensor) may have different sensitivities to ACA and therefore 
provide different estimates of / ACA and (DRE) ACA . For example, one sensor may only 
20 retrieve ACA for optically thick clouds. This sensor would retrieve a larger (DRE) ACA , 
but a smaller f ACA , in comparison with another sensor that is able to retrieve ACA for 
all clouds. Therefore, when comparing the ACA or all-sky DRE estimated based on 
different instruments or algorithms, it is important to compare both f ACA and (DRE) ACA 
terms in Eq. (4). 

25 Acknowledgements. This research is supported by NASA via grant NNX14AB21G. The au- 
thors are grateful to NASA for providing MODIS and CALIOP products used in this study. 
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Table 1 . Quality control metrics used for screening the CALIOP aerosol layer product. 


Criterion 


CAD„score 

HorizontaLaveraging 
Exti nction„QC J532 


< -30 

< 20 km 
0 or 1 


Feature_Optical_Depth_Uncertainty_532 < -99.5 
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Table 2. Regional and seasonal mean values of instantaneous DRE and RFE based on the 
pixel-level computation and the new method. 



DRE [Wm“ 2 ] 

RFE [Wm _2 AOD“ 1 ] 


Bias adjusted 

Bias adjusted 


(unadjusted) 

(unadjusted) 

Pixel computation 

6.63 (5.92) 

55.97 (50.34) 

New Method 

6.39 (5.77) 

53.77 (50.22) 
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Fig- 1- A schematic example to illustrate how CALIOP aerosol layer height information is used 
in our method to determine the population of liquid-phase clouds below the aerosol layer in the 
MODIS COD-CTP joint histogram. 
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Cloud Optical Depth 


Fig. 2 . A schematic illustration of our fast scheme to correct the COD retrieval bias in the 
MODIS cloud product due to overlying aerosol contamination. 
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Fig. 3. (a) Seasonal mean (August/September 2007-2011) instantaneous TOA DRE of above- 
cloud smoke and polluted dust based on the pixel-level computations from Meyer et al. (2013); 
(b) seasonal mean instantaneous TOA aerosol RFE (i.e., DRE per AOD) from Meyer et al. 
(2013); (c) same as (a), but based on the new method; (d) same as (b), but based on the new 
method. 
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Fig. 4. Meridional mean DRE and RFE for the region based on the results in Fig. 3. Lines with 
cross symbol correspond to pixel computations from Meyer et al. (2013). Lines with square 
symbol correspond to results based on the new method. 
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Fig. 5. (a) Annual mean AOD (at 550 nm) of above-cloud light-absorbing aerosols (i.e., smoke 
and polluted dust) derived from 4yr (2007 ~ 2010) of the CALIOP 5 km aerosol and cloud layer 
products, (b) Annual mean below-aerosol COD derived from the MODIS daily level-3 COD- 
CTP joint histogram, (c) Annual mean instantaneous TOA DRE of above-cloud light-absorbing 
aerosols derived using the new method. 
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